Search for: All records

Hydrodynamical interaction in circumbinary discs (CBDs) plays a crucial role in various astrophysical systems, ranging from young stellar binaries to supermassive black hole binaries in galactic centres. Most previous simulations of binary-disc systems have adopted locally isothermal equation of state. In this study, we use the grid-based code Athena++ to conduct a suite of two-dimensional viscous hydrodynamical simulations of circumbinary accretion on a Cartesian grid, resolving the central cavity of the binary. The gas thermodynamics is treated by thermal relaxation towards an equilibrium temperature (based on the constant − β cooling ansatz, where β is the cooling time in units of the local Keplerian time). Focusing on equal mass, circular binaries in CBDs with (equilibrium) disc aspect ratio H/R = 0.1, we find that the cooling of the disc gas significantly influences the binary orbital evolution, accretion variability, and CBD morphology, and the effect depends sensitively on the disc viscosity prescriptions. When adopting a constant kinematic viscosity, a finite cooling time (β ≳ 0.1) leads to a binary inspiral as opposed to an outspiral and the CBD cavity becomes more symmetric. When adopting a dynamically varying α-viscosity, binary inspiral only occurs within a narrow range of cooling time (corresponding to β around 0.5).

Hydrodynamical interactions between binaries and circumbinary disks (CBDs) play an important role in a variety of astrophysical systems, from young stellar binaries to supermassive black hole binaries. Previous simulations of CBDs have mostly employed locally isothermal equations of state. We carry out 2D viscous hydrodynamic simulations of CBDs around equal-mass, circular binaries, treating the gas thermodynamics by thermal relaxation toward equilibrium temperature (the constant-βcooling ansatz, whereβis the cooling time in units of the local Keplerian time). As an initial study, we use the grid-based codeAthena++on a polar grid, covering an extended disk outside the binary co-orbital region. We find that with a longer cooling time, the accretion variability is gradually suppressed, and the morphology of the CBD becomes more symmetric. The disk also shows evidence of hysteresis behavior depending on the initial conditions. Gas cooling also affects the rate of angular momentum transfer between the binary and the CBD, where given our adopted disk thickness and viscosity (H/r∼ 0.1 andα∼ 0.1), the binary orbit expands while undergoing accretion for mostβvalues between 0 and 4.0 except over a narrow range of intermediateβvalues. The validity of using a polar grid excising the central domain is also discussed.